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We propose a model of directed lines where the average direction has the nature of 
a cholesteric liquid crystal. This model, for instance, would describe the liquid of screw 
dislocations in the twist-grain-boundary (TGB) phase of liquid crystals. We show that 
the presence of lines does not alter the long wavelength elasticity of a cholesteric and, 
therefore, does not stabilize Landau- Peierls instability of the cholesteric phase. We discuss 
other possible mechanisms for stabilizing the twist-grain-boundary phase. 
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1. Introduction and Summary 

Directed line liquids have received attention recently in a variety of interesting phys- 
ical situations. Composed of lines which are compelled on average to have tangent vec- 
tors parallel to an imposed axis, they are central to the long- wavelength physics of high- 
temperature superconductors in an external magnetic field, both as entangled flux liquids 
[1] and vortex glass states [2] . They are also of interest in nematic polymers [3] , electrorhe- 
ological fluids and ferrofluids [4]. In this paper, we will investigate a model for chiral line 
liquids in which tangent vectors of the lines are constrained to rotate in a helical fashion 
along some fixed pitch axis. 

The entangled flux liquid in a superconductor is produced when thermal fluctuations 
cause the regular Abrikosov vortex lattice to melt. The close analogy between the smectic- 
A phase in liquid crystals and the superconducting phase in metals and that between 
the nematic-to-smectic-A and normal-metal-to-superconductor transition was pointed 
out some years ago by de Gennes [5]. This analogy leads to the theoretical expectation 
that a liquid crystal analog of the Abrikosov vortex phase should exist with dislocation 
lines replacing vortex lines. Recent experiments [6] confirm that such a phase does exist 
and that its properties are identical to those of the theoretically predicted TGB (twist- 
grain-boundary) phase [7] consisting of periodically repeated twist-grain-boundaries, each 
composed of periodically spaced screw dislocations whose axes rotate in a helical fashion 
from one boundary to the next. The TGB phase can presumably melt to a chiral line 
liquid just as the Abrikosov phase melts to a directed line liquid. The chiral line liquid we 
study in this paper is intended to model the melted TGB phase. 

The entangled flux liquid phase in a superconductor is in fact the normal-metal phase, 
but with enhanced viscosities. Polymer nematics, on the other hand differ from normal 
nematics in that their splay elastic constant K\ is infinite if all polymers extend from 
one end of the sample to the other [8]. Non-linearities constrained by full rotational 
invariance lead [9] to a renormalized theory in which elastic constants take on a momentum 
dependence. In particular, the twist and bend Frank constants diverge at small momentum 
while the splay modulus vanishes. It is thus natural to ask whether there are differences 
between the chiral line liquid and a cholesteric liquid crystal. Our calculations show that 
the long-wavelength properties of the chiral line liquid are identical to those of a cholesteric 
liquid crystal. A cholesteric is a one-dimensional layered structure with a long-wavelength 
Landau- Peierls elasticity and associated destruction of long-range periodic order [10]. The 
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TGB phase is also a one dimensional-layered structure, and though its elasticity is more 
complex than that of a cholesteric [11], its long-range order is also fluctuation destroyed. 
A chiral line liquid has the same long-wavelength Landau-Peierls elasticity as a cholesteric 
but with elastic constant renormalized by the presence of dislocations. The splay elastic 
constant K\, which does not contribute to the Landau-Peierls energy, however, diverges 
as the square of an inverse wavenumber as it does in a polymer nematic [12]. 

In this paper we add explicit degrees of freedom to describe the dislocation line config- 
urations. Though the energy of a smectic-A screw dislocation tilted a small angle a away 
from the director goes as a 2 log(l/a), we will assume that we are at sufficiently long scales 
that the tipping energy becomes analytic [13]. In section II we derive a Landau theory 
for a melted, twisted line liquid, generalizing work on directed line liquids. In section III 
we derive the effective long-wavelength theory and find expressions for the new, effective 
Frank constants of the cholesteric. Finally in section IV we comment on a possible way 
for the non-harmonicity of the tipping energy to stabilize the Landau-Peierls instability of 
the twist-grain-boundary phase. 



2. Derivation of Model 

Our system is a collection of lines described by their locations Rj(s) and parameterized 
by arclength s embedded in and interacting with a cholesteric liquid crystal with a unit 
director field n(r). We model the lines by a local density field p and a local tangent-density 
vector m. In terms of the paths of the lines Ri(s), 

PW = E / ds5 d (R l (s)-r) 7 (2.1) 



and 



m (r) = j ds ^ir 6d{ ^^ - r )- ( 2 - 2 ) 



If our lines do not terminate, then m must be divergenceless: 



da ^(f). w d (R . (s) _ r) 
as 



(2.3) 



= -E/ dsj-6 d (K l (s)-r)=0. 

i 

We will not consider the case of finite lines because dislocations cannot end within a sample 
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In the TGB phase, screw dislocations align parallel on average along the local nematic 
director n. We assume this remains true in the line liquid phase and propose the following 
free energy to describe the chiral line liquid: 

F = J rf 3 r|(m-p n) 2 + F n [n], (2.4) 

where 

F n [n]= l - j d 3 x [^(V-n) 2 + K 2 (n-[Vxn] - k ) 2 + K 3 (nx[Vxn]) 2 j . (2.5) 

is the usual Frank free energy for a chiral liquid crystal and the constraint V-m = is 
understood The free energy favors (m) = p n. We could have included an anisotropic 
tensor coupling of the form 

F = ^Jd 3 x (m M - pon^I^imv - p Q n v ). (2.6) 

By rotational invariance = AS^ U + Bn^,n v . If we keep only terms up to quadratic order 
in the fields, we are left only with the isotropic tensor 8^. In principle, the self-coupling 
of the vector can be different along the pitch axis of the twist-grain-boundary phase than 
perpendicular to it, or in other words, if the pitch points along the x-axis, I xx ^ I yy = I zz . 
We might also include interactions of the defects with themselves involving m by itself. 
We believe, however, that in the long wavelength limit, this model captures the essential 
physics. 

The ground state configuration of the liquid crystal has a spontaneously broken sym- 
metry. We choose the pitch axis to lie along the x-axis and the equilibrium director, 

n (r) = [0, sin(/c x), cos(/co^)], (2.7) 

to minimize the Frank Free energy. We introduce director fluctuations via n = n + 8n. 
If we only work to quadratic order in fields, then because n is a unit vector, it is sufficient 
to consider only 8n with <5n-n = 0. Since (m) is parallel to n in equilibrium, we set 

m = pn + t. (2.8) 

Fluctuations of m along n are described by p and fluctuations perpendicular to n de- 
scribed by t. We, therefore, require t-no = 0. Fluctuations in the equilibrium average of 
m will point along n . Thus (t) = 0, and (m) = (p)n . The constraint V-m = is now 

n -V/9 + V-t = 0. (2.9) 
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and our harmonic theory becomes 



F = |rf 3 r||(t-Mn) 2 + |^ 2 }+F n [n], 



(2-10) 



subject to no-V/O + V-t = 0. We have set p = po + 5 p. and 
F Sn [5n] =F n [n + Sn] 



= 1/^ 



K 1 (V-5n) 2 + K 2 (n -[Vxdn]) 2 + K 3 (k 8nxn + n x[Vxon] )' 



(2.11) 



We note that if ko = then n = z and this theory reduces to the model derived for 
directed line liquids in [12]. 

In order to study the fluctuations in this system we, seek a parameterization of m 
and p satisfying the constraint (2.9) and the constraint n t = 0. The first constraint is 
satisfied by setting 



P = Po- PoV-u 

t = p (n -V)u - /9 (u-V)n , 



(2.12a) 
(2.126) 



where u is any vector. This equation is invariant under the transformation u — > u' = 
u + VxA. This comes about because the constraint (2.9) only involves the longitudinal 
components of t. There is, therefore, a gauge-like symmetry of the transverse components. 
The constraint t-n = implies 

n -V(n -u) = 0. (2.13) 

We can solve this constraint on u by choosing no-u = 0. This is tantamount to a choice 
of gauge in the sense that if we are given a u with no-u ^ we can, without making u 
trivial, choose a A such that n • (u + VxA) = 0. Note that any vector u(r) that depends 
only on x, the component of r along the pitch axis, satisfies the constraint (2.13). Thus 
the free energy expressed in term of u will depend only on gradients of u in the yz-plane. 



3. Rotational Symmetries and Long-wavelength Elasticity 

The free energy for a chiral line liquid, like that for a cholesteric, is invariant with 
respect to uniform translations and rotations. It is precisely these invariances that lead 
to the Landau-Peierls elastic free energy and the destruction of long-range order in a one- 
dimensionally modulated structure. In order to obtain the effective long-wavelength free 
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energy, it is useful to understand how uniform translations and rotations are described by 
the variables u and 5n describing deviations from the equilibrium ground state. Following 
the analysis in [14], we write our fields u and Sn as 

u= [cf),9cos(kox),—9sm(kox)], (3.1a) 
5n= [f,gcos(k x),-gsm(k x)], (3.16) 

thus enforcing both 5n ■ n = and u • n = 0. Because the ground state is not transla- 
tionally invariant, it is convenient to decompose our fields into a sum over Brillouin zones. 
We write for any field X, 

X(r) = J2Xn(r)e ik ° nx , (3.2) 

n 

where X n (r) only has fourier modes with x-components k x G (— 4^-, 

We now consider a rigid rotation described by the vector O. Under this rotation n 
becomes: 

iin =rio + [Cly cos(fco^) — Q z sin(/co^), — O x cos(kox), O x sm(kox)] 

(3.3) 

— (Q> y z — 2 y) [0, ko cos(k x), —ko sin(/c x)] . 

We can now compare (3.1 )and (3.3) to determine what fields / and g are generated by 
uniform rigid rotations. We find 

/ = Vt y cos(Zco^) — Q z sin(/coic) (3.4a) 
g = -Q x - k (Q y z - Q z y) . (3.46) 

If we write t in terms of 4> and 9 and match the same rigid rotation, we find 

Vt y = d z (f) (3.5a) 
-tt z = d y (f) (3.56) 
— fi x — ko(Q y z — 0, z y) = —k (p + sm(kox)d y 9 + cos(kox)d z 9 (3.5c) 

Then writing (3.4) and (3.5) in terms of the decomposition (3.2) gives (with 2X + = 
X, + X_, and 2X_ = -i{X 1 - X_,)) 

f+ = M y (3.6a) 

/_ = Q z (3.66) 

go = -Q x - k {SlyZ - Q z y) (3.6c) 

d y 9_ - d z 9 + + ko4>o = + k (ft y z - Q z y) (3.6d) 

d y 4>o = (3.6e) 

d z ct>o = n y (3.6/) 



d x < 



The free energy cannot depend on the uniform rotation angle O. It will, however, depend 
on linear combinations of the fields 0, 6>, /, and g that do not depend on f2. We can 
construct such combinations with the aid of (3.5) . Note that ft only depends on 9 via 
the combination d y 9- — d z 9 + . The combination d y 9- + d z 9 + does not generate a rotation, 
and the free energy can depend on it. Finally, we note that t will be nonzero if d x cf)Q is 
nonzero. Thus, we expect on quite general grounds that F will have the form 

a [(d y g - 2k f_) 2 + (d z g + 2k f + ) 2 } 
+ b[(d y <p + 2f-) 2 + (d z( Po-y + ) 2 ] 
+ c(ko<po + d y 9„-d z 9 + +g ) 2 }. (3.7) 

+ d(d y 9 + + d z 9_+d x <j> ) 2 
+ e [(d y 9 + ) 2 + (d z 9 + ) 2 + (d y 9.) 2 + (d z 9 + ) 2 ] J 

This form is dictated by the underlying rotational invariance. In terms of the original 
constants in (2.4) and (2.5), a = (Kx + K 3 )/2, and 2b = c = d = 2e = Ap\. The values 
of 6, c, d and e are related by the underlying rotational invariance discussed in section II. 
There are other terms as well, in addition to higher powers and derivatives of the invariant 
combinations. We also note that because p = —p^V • u, <p n for n ^ will be massive. 
Likewise, because of the energy of a splay configuration f n will be massive for n ^ 0. 

If we only consider the cholesteric, then po = 0. In this case we are only left with 
the term proportional to a in (3.7). By integrating out f± we see that d z g and d y g 
cannot appear quadratically. Thus they will first appear quartically, and we find the 
classic Landau-Peierls instability. Including the lines, we see that 9 appears without x- 
derivatives and always in combination with other fields. Upon integrating out 9 both the 
terms proportional to c and d disappear. 

The preceding analysis showed that we must keep all the modes in the first Brillouin 
zone as well as / + and /_. Returning to our proposed model (2.10), we integrate out the 
modes in the second zone, as well as 4>o- We find, to leading order in all derivatives, an 
effective free energy for go : 

'Apl 



d A x 



3k 



(d x go) 2 + ^^[(d 2 y + d!)g \ 



(3.8) 



Thus the fluctuation g suffers from the Landau-Peierls instability. The penetration depth 
of the smectic, A, is 

o 3(K 2 + 3K 3 ) 



8k 2 (3K 2 + Ap 2 ) 7 
6 



(3.9) 



and thus the effect of the defects is to decrease the penetration depth. Though there is 
no long range order, A sets the scale over which it is possible to see the dephasing of the 
ground state. Consider two defects in the same plane of constant phase. Their mutual 
repulsion forces them apart, and as they move they drag the plane with them. Thus the 
effect of the repulsion should be to make the system "less" ordered. 
Additionally, the free energy for the remaining director mode fo is 

K 2 + K 3 



K 1 (d x f y + 



[(d y f ) 2 + (d z f Q ) 2 } + (Ap 2 + K 3 k 2 )f 2 



(3.10) 



and so the fo fluctuations become massive. In analogy with results on polymer nematics 
[12] we could interpret part of this as a defect density dependent divergent shift in K\, 
namely 



K[(a) = K 1 + 



Q 2 



(3.11) 



4. Stabilization of the Twist-Grain-Boundary Phase 

Because the Landau- Peierls instability in the twist-grain-boundary phase is due to the 
underlying rotational invariance, it is hard to imagine how any long-wavelength description 
of the phase could have long range order. As mentioned above, the energy of a defect tilted 
at an angle a with respect to the layer normal is proportional to a 2 log(l/a). While it 
is true that upon integrating out short-distance modes, the effective energy will become 
analytic in a, the coarse-graining procedure will be cutoff by length scales inherent to 
the system. In particular, as we coarse-grain, the distance between the screw dislocations 
decreases. When the distance between them becomes on the order of the penetration 
depth, we must stop. If the screw dislocations are close enough (i.e. ko is sufficiently 
large), then the defect energy will remain non-analytic. The argument in the previous 
section would break down, because now we could have a term such as — log((/)t(q) -t(— q) 
in the free energy in addition to those in (2.4). Though again, we would expect that t 
would always appear in combination with 5n, it is hard to see how non-analytic momentum 
dependence would change things. Since the energy cost of a bend in the direction of the 
cholesteric pitch will pick up this logarithmic energy in the twist-grain-boundary phase, 
we might consider, as a toy model, a slightly more rigid cholesteric described by a phase 
field u with free energy 

F= \j (03 W (- c l)[ S1 °gl 1 /^x|^ + ^ + ? 2 2 ) 2 ]«(q). (4.1) 

7 



In this model, there is still a Landau-Peierls instability, in that 




f d 3 q 1 



J (2tt)3 -B\og\aq x \q*+K(q2 + q2)2 
f dq x q±dq± 1 



J (2tt)2 -B\og\aq x \ql + Kqi 



(4.2) 



1 f 1/a dq. 



32tv^/KB Ji/l q x V- log aq, 



y/log(L/a) 



and so there will still be a divergence with system size, but it will scale as the square-root 
of the logarithm instead of the logarithm itself. 

Unless there are sufficiently long-ranged interactions it appears that the underlying 
rotational invariance will always destabilize the twist-grain-boundary phase. 
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